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Abstract 

We propose a simple model for Warm Dark Matter (WDM) in which two fermions 
are added to the Standard Model: (quasi-) stable "keVins" (keV inert fermions) 
which account for WDM and their unstable brothers, the "GeVins" (GeV inert 
fermions), both of which carry zero electric charge and lepton number, and are 
(approximately) "inert" , in the sense that their only interactions are via suppressed 
couplings to the Z. We consider scenarios in which stable keVins are thermally 
produced and their abundance is subsequently diluted by entropy production from 
the decays of the heavier unstable GeVins. This mechanism could be implemented 
in a wide variety of models, including Eq inspired super symmetric models or models 
involving sterile neutrinos. 
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1 Introduction 



One of the biggest questions in contemporary astroparticle physics is about the nature of 
Dark Matter (DM), which makes up more than 20% of the energy content of the current 
Universe pp. When considering particle physics motivated candidates for a DM particle, 
one distinguishes between cold (CDM) and hot Dark Matter (HDM), depending on if the 
freeze-out of the corresponding particles happened when they were already non-relativistic 
or when they still had relativistic velocities. Since a possible DM candidate particle should 
be electrically neutral, not strongly interacting, and (quasi-) stable, the only possibility 
within the Standard Model (SM) of particle physics would be a light neutrino. Due to 
strong experimental constraints (from single- f3 decay experiments [21 [3], double-/? decay 
experiments [H [5], and cosmological bounds [1]), the neutrino mass should be less than 
about 1 eV, making it a perfect HDM candidate. Unfortunately, HDM is excluded as 
the dominant DM component due to problems with structure fomation 0, [7J. On the 
other hand, many theories beyond the SM predict so-called WIMPs (weakly interacting 
massive particles), which would be very suitable CDM candidates. 

However, recently the intermediate case of warm Dark Matter (WDM) has attracted 
considerable attention, in particular from the structure formation side. After pioneering 
works [H |9] more than ten years ago, several simulations have been performed (see, e.g., 
Refs. [TOl ELI]), in particular to study the small scale structure [121 HSl EH] induced by 
WDM, in order to maybe be able to solve the typical CDM problem of predicting many 
more dwarf satellite galaxies than observed. In addition, model-independent analyses [151 
[To] and surveys such as ALFALFA [T7J seem to point into the direction of WDM. In 
general we shall refer to the candidates responsible for WDM as "keVins" (keV-mass 
inert particles). Just like WIMPs, keVins also share certain common properties, which is 
why it makes sense to identify them as a general class of WDM candidates that can be 
studied both in a general context and in a specific model. 

Probably the leading particle physics candidate for a keVin that has been studied 
most frequently is a sterile neutrino with a mass of a few keV. Frameworks in which such 
a particle can appear are, e.g., the so-called i^MSM [18], the left-right (Li?)-symmetric 
framework [19], the scotogenic framework [201 |2T] . or they could appear as composite 
neutrinos [22]. When studying keV sterile neutrinos as DM candidates, there are two 
main tasks, namely finding a suitable production mechanism to generate the correct DM 
abundance and explaining the very existence of the keV mass scale. Possible production 
mechanisms are, for example, non-resonant [23] and resonant [21] non-thermal production 
or thermal production with subsequent entropy dilution [25] , which all have been applied 
to the keV sterile neutrino case P31 ESI EZl EBl ESI EQl EH E2J - To motivate the existence of 
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a keV scale in the first place, the possibilities include the split seesaw mechanism [33"1 134]. 
flavour symmetries [351 ESI EI] , the Frog gatt- Niels en mechanism [381 ESI SO] , or extensions 
of the seesaw mechanism [U]. See, e.g., Refs. [121 BEl 133] for more general discussions. 
However, apart from keV neutrinos, there could also be other particles playing the role 
of keVins. The cases studied in the literature include the gravitino [351 SHI E2] 5 the KK- 
graviton [46J, singlinos [H], axinos [IS], axion-axino mixed DM [39], or majorons [50l [5Tj . 

In this paper we propose a new and very simple model of WDM based on adding a pair 
of electrically neutral fermions to the Standard Model: a stable keVin, xi> accompanied 
by an unstable partner, the "GeVin" X2, with a mass in the GeV range. Both particles in- 
teract only by suppressed couplings to the SM gauge boson Z. These suppressed couplings 
may arise, e.g., from a small mixing effect in more sophisticated models. We show that 
the correct Dark Matter relic abundance can be obtained by thermal overproduction of 
Xii which is subsequently diluted by the production of additional entropy from the decay 
°f X2-, while satisfying constraints from structure formation and big bang nucleosynthesis. 
Such a mechanism could be applied to a wide variety of models, for example E§ inspired 
supersymmetric models [52] or models involving sterile neutrinos [19]. However, we stress 
that it is also possible to have keVins without GeVins, with the keVin being produced 
non-thermally due to its feeble interactions, which is the traditional mechanism for warm 
Dark Matter from keV sterile neutrinos. 

This paper is organized as follows. In Sec. [21 we first introduce the general framework 
of keVins and perform a detailed semi-numerical analysis of how to generate the correct 
DM abundance using thermal (over-) production of the keVins in the early Universe. 
This is achieved by suppressed couplings of the keVins to the Z-boson. The abundance is 
later on diluted by additional entropy produced in decays of a heavier electrically neutral 
particle, the GeVin \i- We also identify several benchmark points for the five independent 
parameters, which can then be used to check specific models. Exactly that will be done 
in Sec. El where we first present a situation in the framework of the E^SSM in which 
keVins naturally arise. After that, we also give a short discussion on how to extend our 
consideration to the case of keV sterile neutrinos in a Li?-symmetric framework. These 
examples serve as guidelines on how to apply our considerations, but they are certainly 
not the only examples one could find. Finally, in Sec. HI we will summarize our findings 
and conclude. 

2 keVins and GeVins 

In this section, we introduce a simple model of WDM in which two fermions are added 
to the Standard Model: keVins and GeVins, both of which carry zero electric charge and 
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lepton number, and only interact via suppressed couplings to the Z. Later on, in Sec. |3j 
we will apply this mechanism to other examples, e.g. within the i? 6 SSM or to models 
with sterile neutrinos. We want to stress, however, that such examples do not play any 
decisive role, and that other examples could easily be found which also involve the simple 
mechanism discussed in this section. 

2.1 The Z couplings 

The starting point of the model is to add a pair of electrically neutral fermions to the 
Standard Model: a stable keVin, Xi, accompanied by an unstable GeVin, xi- Both these 
fields are assumed to have suppressed couplings to the SM gauge boson Z. This may be 
achieved by having a singlet fermion with small mixing to other fermions, as in the E$ 
models discussed later in the paper. However, here we simply assume suppressed couplings 
of keVins and GeVins to the Z. Note that the electric charges of Xi are strictly zero. In 
addition, we consider the possibility for the Z to also couple in a flavour-changing way to 
Xi and X2, allowing decays such as xi ~ t Xi + Z*, Z* ~ * f f ■> where / is any kinematically 
accessible SM fermion. 

The masses of Xi,2 are denoted as M 12 . What we have in mind is to have masses 
Mi ~ 0(l - 100 keV), while M 2 > 0(1 - 100 GeV) is considerably heavier. Hence, we 
have arrived at a framework where we have one particle with a keV mass that is stabilized 
by some additional quantum number, for example .R-parity in supersymmetric models. 
This particle is mainly a SM singlet fermion, which is why the name keVin (keV-mass 
inert particle) is appropriate. Furthermore, the keVin is the WDM candidate. 

Before we give an estimate of the DM-abundance, which turns out to be easily and 
naturally reproducible in this simple framework, we want to comment on the bounds 
arising from the Z-boson decay width: both our singlets, the keVin xi as we U as its 
heavier brother, the GeVin Xii wm have suppressed couplings to the SM Z-bosons, given 
by the small parameters and 5. We assume the following interaction Lagrangians|j 

• Z-x t -Xi (* = M): 

Cu = gZ^lxil^lsXh (!) 

£12 = geie 2 SZ ll xil' 1 'y5X2 + h.c, (2) 

1 Note that we assume for the x< a pure axial vector coupling to the SM Z-boson, for simplicity. 
Our considerations would, however, not be considerably altered if any type of mixed vector/axial vector 
coupling was assumed, and the reader is invited to repeat our analysis for such a case. 
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where g = 0.653 is the usual SU{2) gauge coupling, and ^ are the suppressed couplings 
Xi to the Z. Note that, if the masses M 12 are smaller than Mz/2, they would contribute 
to the invisible Z-boson decay width, whose value is given by r^ invisiblc = (499.0 ± 
1.5) MeV [53]. Since this value is more or less exactly what is obtained by considering 
the decay of Z into pairs of the three known active neutrino flavours, the corresponding 
decay width into inert fermions, 

T x = T(Z xiXi) + nZ X1X2) + T(Z X2X2) = (e{ + e\ + 2ej 2 ) , (3) 

where ei 2 = €ie 2 S, must be smaller than roughly the uncertainty of the invisible decay 
width, T x < 1.5 MeV. Assuming that a certain et dominates, this will turn into a bound 
of 3 < 0.23. 

This leaves us with only five decisive parameters: M 12 , ei^, and S. 
2.2 Freeze-out regions 

The crucial values to be determined are the freeze-out temperatures of the two particles 
under consideration, xi an d Xi- Therefore, we consider the annihilations of a Xi~Xi P arr 
into either a SM fermion-antifermion pair, or into a pair of VT-bosons, by a Z-boson 
mediator, depending on the temperature and initial state mass. Note that we neglect 
meson final states since they will, in the relevant parameter regions, either be negligible 
or kinematically not accessible. The corresponding spin-averaged matrix elements 
and |.M|^ for annihilation into fermion pairs or VT-bosons, respectively, are given by 

| - M| ' = c^[{s-Miy + MlTlf K^)(*aP2) + (k lP2 )(k 2Pl ) - M?^)] , 

W&r = {s _M2$f M 2 T 2 Hpi ~ ft.) " HPi - Pa) + ipm){hk 2 - Mf )] . (4) 

Here, hi (p^ are the initial (final) state momenta, cw = cos8w is the cosine of the 
Weinberg angle, T z = (2.4952 ± 0.0023) GeV is the Z-boson total decay width, y/s is 
the center-of-mass energy, and S = J2f {A) + Bf) with A f and B f given in terms of SM 
fermion vector and axial vector couplings as Af = gy + g^ and Bf = g v — g^, where 
the sum runs over all fermions that are kinematically accessible. We have plotted £ as 
function of the temperature in Fig. [T], left panel. 

The two limiting cases for which analytical solutions to the Boltzmann equation ex- 
ist are those of relativistic freeze-out (hot thermal relic) and of non-relativist ic freeze-out 
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(cold thermal relic). For both these regions, we have first computed the corresponding ap- 
proximations of the expressions in Eq. (jj]), and have then performed the thermal average. 
In case of relativistic freeze-out, this results in 

where F(T) ~ A ( r^ijzjs ) M(l)> an d the function A(X) is given by 

oo 

A(X) = J Xs^K, . 



(6) 



s=0 



Note that A(l) = 32 T 5 can be calculated analytically, while A (^ jjzm^s ) * s evaluated 
numerically. In the non- relativistic limit, in turn, we obtain^ 



9 4 efX MjT w g*<*c? w M{T 

2^ w ' (s - M\f + Mfr| ana i. 47r • {s _ M | )2 + M |r|- 



(7) 

Furthermore, we need the number densities for the relativistic and non-relativistic cases, 
which are given by 



3 C(3) ^3 A (MiT\ 

n m l. = \9 X —^T T aIld ^non-rel. = 9 X \ J 



3/2 



e -Mi/T 



where g x = 2 is the number of internal degrees of freedom of the fermion x%> C( x ) is 
the Riemann zeta function, and we have neglected the chemical potential in the non- 
relativistic case. 

The freeze-out temperature is then determined by setting the product n(av) equal to 
the Hubble rate H = jj^^9~* where M P = 1.22 ■ 10 19 GeV is the Planck mass 
and g* is the effective number of relativistic degrees of freedom (cf. Fig. [I] right panel). 
Finally, by setting the freeze-out temperature Tpo in the respective case equal to the mass 
Mi, we can obtain the (approximate) range of validity for both limiting cases. The result 



2 Note that, due to having two Majorana fermions in the initial state, the leading term in the non- 
relativistic approximation of the annihilation cross section is the p-wave, which is suppressed by a factor 
v 2 of the velocity. This translates into a factor 6T/M2 in the thermal average. For a detailed discussion 
of the cases where such suppressions are present or not, see Ref. [53] . 

3 There is no need to perform a more sophisticated procedure to determine the freeze-out temperature, 
since we will in any case use some other approximations, in particular when solving the Boltzmann 
equation. 
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Figure 1: The quantities E and g* as functions of the temperature. Thereby, g* includes 
the contributions from the two light fermions that are present in our framework in addition 
to the SM particle content, i.e., the keVin xi an d the GeVin \2i making it different from 
its SM counterpart. Note that the biggest "jump" in g* is at the QCD scale, where 
quarks are bound into nucleons. We have taken the example value of Mqcd = 300 MeV 
for defmiteness (see Ref. [53] for more details on that point). 



is displayed in Fig. [2j In the two limiting regions, it is possible to obtain an analytical 
estimate for the DM number density normalized by the entropy density, Y Xi = rii/s, while 
in the intermediate region one would need to numerically solve the Boltzmann equation to 
obtain Y Xi . Such a numerical solution is, however, beyond the scope of this paper, and we 
will leave that more elaborate investigation for future work. Hence, in our case, we obtain 
an additional "constraint" on the masses. This does not play a role for M\ ~ keV, which 
will always be in the relativistic freeze-out region. The heavier particle X2, in turn, could 
freeze-out relativistically or non-relativistically, depending on the exact combination of 
mass M 2 and coupling suppression e^- 

Note that, for a too strong suppression 6j in the annihilation cross section, the product 
n(av) will always be smaller than H, i.e., the particle never enters thermal equilibrium. 
In that case, thermal freeze-out is not possible, which is marked by the left gray rectangle 
in Fig. For too large couplings, in turn, and a mass Mj < M z /2, there is the bound 
by the Z-boson decay width, marked by the right gray rectangle in the plot. 

4 It may, however, be possible for very small couplings to produce the correct DM-abundance by the so- 
called freeze-in instead. In that case, the corresponding particle is often called FIMP (feebly interacting 
massive particle) |55j . 
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Figure 2: Regions of relativistic and non-relativistic freeze-out. The intermediate region 
is marked by the light gray band. (See text for further explanations.) 

2.3 Analytical estimates 

Before turning to a detailed analysis, it is useful to give some analytical estimates. The 
principle idea is the following: since the keVin xi freezes out while still being relativistic 
and since it has a mass Mi of a few keV, its abundance will be much too large for it to play 
the role of the DM. However, if there is at some point after the freeze-out of Xi a phase 
where the energy density of the Universe is dominated by the heavier fermion Xii t ne ou t- 
of-equilibrium decays of the heavier particle (into relativistic final states, i.e., radiation) 
could produce a certain additional amount of entropy, parametrized by the factor S, since 
the entropy of radiation is much higher than that of matter. This additional entropy 
delays the cooling of the Universe [25] and it dilutes the natural abundance of xi by 
exactly that factor S [191 126], which can correct its value in such a way to meet the 
observed region Jf| 

5 Note that this argumentation does not take into account the additional population of xi by the 
decays X2 — > Xiff ( an d X2 — > XiW~W + for very large M2). This should be a very good approximation: 
even though X2 dominates the energy density of the Universe, ptot — p(X2)i a t the time of its decay, its 
number density n(x2) — p{X2)l^2 is much smaller than the number density of Xi, due to M2 S> M\ 
and the late freeze-out. Thus, the main contribution of this decay to the DM abundance is indeed the 
dilution by entropy production. 
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In that case, the final xi abundance will be given by [19] 



7.61 -10 2 Mi g x , eS 
Xl ~ S lOkeV gUTjoV [ } 

where g x>e s = \g x an d 9*s(Tp ) ~ g*(Tp Q ). Here, Tp Q denotes the freeze-out temperature 
of Xi- The factor S is given by the ratio between the entropy Sf after most of the X2 
decays [i.e., for t > t(x 2 )] and the entropy Si before the decays set in. This ratio is 
derived in detail in Ref. [23], and we will essentially make use of Eq. (23c) therein: 

where we have approximated the weighted average of the effective number of relativis- 
tic degrees of freedom by the actual such number g* at the time under consideration. 
Furthermore, Y" 2i00 is the final abundance of X2 particles as obtained by the Boltzmann 
equation in the limit of small temperatures, and T 2 is the decay width of X2- The decays 
under consideration involve a flavour changing coupling of the Z-boson, so that X2 decays 
to Xi an d a fermion-antifermion pair. |^| Taking the final states to be effectively massless, 
it is easy to derive the decay width: 

r(xa -> xrfl) ~ ^ (11) 

To obtain an easy estimate, we should first note that the natural value of the DM 
abundance, Eq. ([9]) with S = 1, is larger than the range allowed by WMAP-7 p], fl-DMh 2 = 
0.1120 ± 0.0056, by a factor much larger than 1. A large correction factor S is achieved 
most easily if X2 freezes out while still being relativistic, where we have Y 2:O0 = ~jf ; 
since in that case the number density of X2 is not suppressed. Hence, in the region where 
we could potentially reproduce the correct abundance, one can neglect the first term in 
Eq. (TiTJl) [T9] , which allows us to write 



S w 0.76 ^ y * , 12 

2 <?*VTWp 

where we have again approximated the averaged value of g* by an example value in the 



6 Note that, in principle, there exists also a decay mode \i — > 3xi- However, this mode is suppressed 
by yet four more powers of ei, so we neglect it here. In principle X2 can also decay into xi phis a pair of 
VF-bosons via the Z-W-W coupling, however over most of the physical region of interest this decay is 
kinematically forbidden. 



8 



region of interest, since it is slowly varying anyway in most of the parameter space. In 
order to obtain the required decay width r req needed to meet the correct relic abundance, 
we can combine Eqs. (Q and (fT2|) : 

In our case, we have g* = g*(Tp ) (v e , M ,r; e± ; an d Xi are reiativistic)0 and hence 

Mo \ 2 /keV x 5 



r,^4.1.10-GeV^7(^j ^-j . (14) 

There is an important bound on the decay width T 2 to be taken into account: in order 
not to disturb Big Bang Nucleosynthesis (BBN), the lifetime of the Xi decay should be 
small enough, i.e., it should be smaller than something like 1 sec [19] . More precisely, the 
reheating temperature, 

r r -if^V 1/1 v / fWp, (is) 



2 V 45 

due to Xi decays should be larger than about 0.7 MeV [56] (weak BBN bound) to 
4.0 MeV [[57j (strong BBN bound), depending on what is taken into account in the analy- 
sis. If this is the case, then the Universe simply continues to cool down after the entropy 
production has been finished, and will at some point undergo BBN in the usual wayjfl In 
our case, this translates into a minimal value of the decay width of 

r min = v 7 ^ 3 - 5 • 10 ~ 24 GeV > 11 ■ 10 ~ 25 GeV )> ( 16 ) 

for T r = (4.0, 0.7) MeV. One can see from Eqs. (fT4|) and (fl6|) that the "natural" value of 
r req (for Mi = 1 keV and M 2 = 1 GeV) is actually below the bound r min . However, this 
is not a problem as long as we can choose the masses Mi and M 2 such that we are okay 



7 Note that this leads to g* w 10, which is the smallest possible value of g* before the inset of Big Bang 
Nucleosynthesis. However, if we choose a smaller coupling e\ it may be that the keVin xi freezes out 
much earlier, at a higher temperature. If this freeze-out happens before the QCD phase transition, while 
the entropy producing decays of X2 happen after the QCD transition, then the factor [.9*(^fo)/s*] 2 leads 
to a significant boost of the RHS of Eq. (jT3|) , which could, for t\ = 0.001, be as large as about 25. This 
would lead to an increase of the reheating temperature by a factor of about -\/25 = 5 and to a reduction 
of the natural xi abundance by about g*(Tp Q ,ei = 0.001)/g*(Tp O , t\ = 0.2) ~ 5. Even though this will 
increase the lower bound on M\ it could still lead to a better situation in total. To take into account this 
effect, we will always consider that the entropy producing decays happen very late, i.e., when g* « 10. 

8 Note that there is one additional relativistic degree of freedom during BBN, namely the keVin xi- 
However, since the abundance of this particle gets diluted by the entropy production, this is not going to 
have a big influence on BBN. 
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with all bounds. 

Requiring r req > r min , one obtains the relation 



> (9.2,1.6)* (17) 



GeV v ' ' keV 

Successful structure formation (i.e. the correct velocity distribution of the DM particles) 
is probed by the so-called Lyman-a (Ly-a) bound [ID] . As explained in Ref. (TH], we 
have to rescale the bound of 8 keV obtained in Ref. [TU] for non-resonantly produced 

— 1/3 

WDM particles by a factor iS rcq , where iS rcq is the entropy dilution factor required to 
meet the correct DM abundance. The corresponding lower bound on M\ is displayed in 
Fig. El where we have used the two example values €\ = 0.001 and e± = 0.2, illustrating 
the robustness of the bound. Since the latter value is very close to the upper bound on e\ 
from the Z-boson decay width, we can conclude that the lowest possible value for M\ is 
around 1.5 keV, which is practically identical to the bound obtained in Ref. [19] . Hence, 
an absolute lower bound on M 2 can be obtained: 

M 2 > (13.8 GeV, 2.4 GeV). (18) 

Glancing at Fig. [2], one can see that the more restrictive (strong) bound pushes the allowed 
region of M 2 and e 2 into a very small region in the upper left corner of the relativistic 
freeze-out region. This does, however, not mean that the whole mechanism does not 
work, but it means that we would need to do a full numerical analysis of the Boltzmann 
equation to treat the intermediate freeze-out region properly. On the other hand, the less 
restrictive (weak) bound leaves us with a much larger region that we can calculate, where 
e 2 can be as large as approximately 5 • 10~ 3 . 

To give an estimate for our parameters, we can use the low-energy approximation of 
Eq. pip as estimate for the full decay width, 

r 2 ~ iS™* 4 B 1.1 . 10-" GeV e f-^-V. (19) 
Translating this into a value for £ = eie 2 5, one obtains 

^2.0.10-W^Xf 2 ^. (20) 

Hence, e.g. for M\ = 2 keV and M 2 = 14 GeV, one obtains £ = 3.6- 10 -9 . Mi is consistent 
with the Ly-a bound for €\ = 0.2, and the value of M 2 restricts e 2 to be something like 
8 ■ 10" 5 , in which case we obtain 5 = 2.2 • 10 -4 . We observe that, because the allowed 
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Figure 3: The Ly-a lower bound on M 1; as function of the natural abundance for certain 
parameter combinations. A working (green) and an excluded (red) parameter combination 
are indicated. 

range of e<i is quite restricted, we need a strong hierarchy between the two suppressions, 
6i < 6 2 , for relatively large values of e\. 

The estimate for the situation with relativistic freeze-out of \2 is also depicted in Fig. HI 
Using the analytical estimates from above, one can already see that the allowed regions for 
entropy production factors of 100 or more are not very big. This tension will even become 
stronger when taking a closer look: in Fig. HI the region left of the vertical line is where a 
relativistic freeze-out of X2 could be possible. However, if it does indeed happens will also 
depend on the value of the coupling constant €2 itself, and we will see that this will in the 
end make it hard to conform with the hard BBN bound on the reheating temperature. 
However, this bound can still be avoided in some circumstances [57], which depends on 
the details of the reheating process. A dedicated investigation of this aspect is beyond 
the scope of this paper, but we would like to suggest it as possible further study. 

2.4 The numerical analysis 

We will now perform a more detailed analysis of our model. For certain scenarios, i.e. 
certain values for the parameters (e\,5), we have first calculated the freeze-out tempera- 
tures, then determined the natural Xi an d Xi abundances, calculated the possible amount 
of entropy production (for certain combinations of M2 and 62), and finally obtained the 
diluted xi abundance for a certain value of M±. Note that, although Q Xl h 2 oc Mi [cf. 
Eq. (j2J)], Mi is still essentially unconstrained by the other parameter values, as indicated 
in Fig. [2j Hence, as long as we can avoid the Ly-a bound, cf. Fig. [3J we can more or 
less select between a range of different choices for M\. This will of course affect the 
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M 2 [GeV] 

Figure 4: Approximate regions where potentially much entropy could be produced, for 
the case of relativistic freeze-out of \i- Note, however, that the vertical line only denotes 
the value of M 2 below which relativistic freeze-out could be possible at all, while the exact 
bound also depends on the value of e 2 (cf. Fig. [2]). 

final Xi abundance, but it will not affect the amount of entropy production, so that in 
one and the same scenario, i.e. for the same combination (e\,8), one might be able to 
hit the observed value of the DM abundance for relativistic freeze-out of X2 (small e 2 to 
have early freeze-out, but larger M 2 to have a fast enough X2 decay) with larger Mi, or 
for non-relativistic freeze-out of xi (large e 2 to have late freeze-out, but smaller M 2 to 
have a slow enough xi decay) with smaller M\. The lesson to learn is that, although 
one might naively favour relativistic freeze-out of X2 to have a high enough abundance, it 
turns out that with non-relativistic freeze-out it is much easier to avoid the BBN bound, 
so depending on the actual choice of parameters one or the other situation could be of 
advantage. 

Let us first have a look at the correponding reheating temperatures, which are depicted 
in Fig. Here, we have calculated Eq. (TT5|) for different combinations of the important 
parameters. The colour code is that red lines signal relativistic (hot) freeze-out of the 
GeVin X2, while blue lines stand for non-relativistic (cold) freeze-out of X2- The interme- 
diate freeze-out region (cf. light gray band in Fig. [2]) has not been calculated by us, due to 
the lack of an approximate solution of the Boltzmann equation, but it is indicated in the 
plots by the thin dashed black lines. However, these lines simply interpolate between the 
hot and cold freeze-out regions, and they can only serve as guide for the eye but do not 
signal the physical reheating temperatures in that region. As we can see from the plots, it 
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Figure 5: Reheating temperatures after x 2 -decay. The gray patches mark the BBN bound 
(dashed: hard bound, T r > 4.0 MeV; solid: weak bound, T r > 0.7 MeV). (See text for 
further explanations.) 



is extremely difficult to have a high enough reheating temperature for relativistic freeze- 
out of \2i while this is relatively easily achieved for non- relativistic freeze-out. However, 
this does not yet give information about how much entropy can actually be produced in 
the X2 decays, which is what we will have a look at next. 

The allowed values for the entropy dilution factor, for the same scenarios as in Fig. 
are depicted in Fig. EJ We have plotted the amount of entropy produced in the \2 decays 
[i.e., the factor S in Eq. (ITU]) ] for different values of M 2 and e 2 , where T = 4.05 denotes the 
thermalization limit leading to the smallest possible value of e 2 (cf. Fig. [2]). As already ex- 
plained earlier, we have only calculated the cases of relativistic (red) and non-relativistic 
(blue) freeze-out of X2, whereas the intermediate region is only indicated by the dashed 
light gray lines (which are again only a guide for the eye interpolating between the two 
limiting cases). Nevertheless, due to the continuous nature of the functions involved, for 
cases where there is some allowed region for the entropy production by relativistically 
frozen-out X2-particles, there must also be some significant amount of entropy production 
in the intermediate region, although the true behaviour may not exactly follow our inter- 
polation. Since we, however, wanted to give a proof of principle first, we leave the more 
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detailed numerical calculation for further studies. 

As already explained, the amount of entropy that can be produced is limited by 
the bound from BBN. In Fig. [6j dark coloured dashed lines denote the regions that are 
consistent with the hard BBN bound (T r > 4.0 MeV), while the solid lines denote the 
regions consistent with the weak BBN bound (T r > 0.7 MeV) only. The lines drawn in 
light red and light blue are excluded by BBN. As a general tendency, we see that it is 
easier to produce enough entropy for the case of relativistic freeze-out of \2 (dark red 
lines), which is exactly contrary to the situation for the reheating temperature in Fig. |5j 
This can be understood by the suppression of the number density for the case of cold 
relics. However, a large mass M 2 can partially compensate for that, as the abundance 
itself is the decisive quantity in Eq. ({TO]) , which is why for some choices of parameters 
there can also be a considerable entropy production for cold freeze-out of the GeVin (cf. 
right panels of Fig. [6]). However, in general the strong BBN bound is a very restrictive 
limit that does not allow for much entropy production, as it forces the decay width T 2 to 
be large. On the other hand, if one only requires the weak BBN bound to be fulfilled, it 
is indeed possible to produced entropy correction factors up to about S ~ 100. 

In principle, it would be possible to produce considerably more entropy by introducing 
more "generations" of GeVins: if there were N particles with a mass similar to M 2 and an 
interaction strength similar to the one of the GeVin, then the entropy correction factor 
would be multiplied by exactly that number, S — >■ N ■ S (if S 3> 1). This would not 
influence the reheating temperature, as it only depends on the value of T 2 , cf. Eq. (fT5l) . 
and it would even strengthen our assumption of the GeVin temporarily dominating the 
energy density of the Universe. Hence, introducing more types of GeVins could be a 
natural way to increase the amount of entropy produced considerably. In particular, one 
could achieve the correct DM abundance while being in accordance even with the hard 
BBN bound. However, although this is a tempting possibility, we will stick to the minimal 
situation in order to present an analysis of the simplest framework. 

The results for the diluted abundance, Eq. (ITU1) . are plotted in Fig. Again, we have 
used the color code of red lines signaling relativistic freeze-out of \ 2 while blue lines signal 
non- relativistic freeze-out of \ 2 - We further distinguish between parts of the parameter 
space that are okay with the strong BBN bound (dark solid lines), parts that are only 
consistent with the weak bound only (dark dashed lines), and regions that are inconsistent 
with both of them (light colored lines). The exception to this is the intermediate freeze- 
out region, which we did not investigate (again, the corresponding gray dashed lines 
are merely interpolations between the hot and cold regions and serve as guide for the 
eye only): in case that already the non- relativistic freeze-out region is inconsistent with 
both BBN bounds, we can be sure that also the intermediate region will not cure that. 
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Figure 6: Entropy production for different values of (M2, 62) in a given scenario, for the 
cases of relativistic (red) and non-relativistic (blue) freeze-out of X2 (the interpolation 
for the intermediate freeze-out region is indicated by the dashed gray line). Dashed 
dark lines are consistent with the weak BBN bound, while solid lines are even consistent 
with the strong bound. Everything printed in light colours is excluded by BBN, i.e., the 
corresponding regions are inconsistent even with the weak bound. 



However, for the other cases we cannot be sure if they agree with the BBN bound, or 
not. To investigate those regions in detail we refer to a dedicated numerical study of the 
Boltzmann equation. 

Apart from that, we can clearly see from the plots what the effect of the entropy 
production is: in case that the factor S is significantly larger than 1, the abundance 
of the keVin xi is diluted by the entropy production, and hence decreased compare to 
its natural value represented by the horizontal black dashed line in Fig. [7J Depending 
on the exact values of M2 and 62, this dilution can be strong enough for the actual xi 
abundance to meet the 3cr range allowed by WMAP-7 p] (green horizontal band in the 
plots). Note that, although one would naturally expect the relativistic freeze-out of X2 
to be the favoured situation (as in the left panels of Fig. [JJ for values of €2 close to the 
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Figure 7: Diluted Xi abundance for different values of (M 2 , e 2 ) in a given scenario, for 
the cases of relativistic (red) and non-relativistic (blue) freeze-out of \2- The natural 
(undiluted) abundance is marked by the horizontal black dashed line, and in general solid 
lines mark consistency with both BBN bounds, dashed lines mark consistency with the 
weak BBN bound only. Light colours signal that the corresponding parameter region 
is inconsistent with BBN. Gray dashed lines correspond to the intermediate freeze-out 
region. (See text for further explanations.) 



thermalization limit), a suppression of the \2 decay by a small parameter 5 can be easily 
cured by a slightly larger value of the mass M 2 , thereby leading to a more advantageous 
relation between the decay width T 2 and the energy density Y 2tO0 M2 in Eq. fflOl) . This 
results into efficient entropy production even for non-relativistic freeze-out of \i (cf- lower 
right panel of Fig. [7]). Note that a further dilution of the natural abundance of the keVins 
Xi can come from a large g* in Eq. (Q: for very early freeze-out of %i, there are more 
relativistic degrees of freedom available in the Universe, and all of them will be produced. 
By this, the energy density present is distributed among more relativistic species, leaving 
less left for xi- Finally, as explained before, we have chosen the values of Mi in order to 
be consistent with the Ly-a bound, but since they are otherwise essentially unconstrained 
one could also choose different values. 
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Summing up, one has many possibilities to indeed hit the correct abundance in our 
setting, which leads us to the conclusion that it would without any doubt be worth to do 
a dedicated numerical study in order to investigate different scenarios in greater detail. 
We also plan to perform such a study ourselves in the foreseeable future. 

2.5 Benchmark points 

Let us look in more detail into the promising parameter regions. We have seen from Fig. [7] 
that in general smaller values of e\ seem to be preferred. This was to be expected, since 
smaller e\ leads to an earlier freeze-out of Xi an d hence to a larger value of g*, which 
decreases the Xi abundance according to Eq. Q. Since it is non-trivial to produce large 
amounts of entropy dilution, i.e. S > 1, a naturally small abundance of xi is desired. 
Furthermore, it is of advantage to "tune" the decay width T 2 separately by adjusting the 
parameter 5, in order to keep the reheating temperature high enough while staying close 
to the lower bound in order to produce enough entropy. Still, from the lower panels of 
Fig. we can see that it is possible to hit the correct abundance for both, hot and cold 
X21 an d in accordance with the weak BBN bound. 

This can be seen more clearly in Fig. [HI where the possible freeze-out temperatures 
(upper panels) and the possible amounts of entropy dilution (lower panels) are depicted, 
for both the HOT (left panels) and the COLD (right panels) scenarios, defined in the 
figure, which correspond to relativistic and non-relativistic freeze-out of X2, respectively. 
This time, we have decided for a linear scale in order to illustrate that the region in 
parameter space leading to significant entropy production is not infinitely small. Note 
that in particular for the HOT scenario, it might very well be possible to find patches in the 
parameter space corresponding to the intermediate (warm) freeze-out region of X2-, which 
could still be in agreement with all bounds and lead to a significant entropy production. 
Overall one observes that it is difficult to simultaneously produce large amounts of entropy 
while still being in agreement with the hard BBN bound, in accordance with our findings 
from Sec. 12.31 However, the weak BBN bound can easily be fulfilled, and in any case a 
very detailed derivation of this bound, preferably within the framework of our model, is 
currently not available. This offers a great possibility to falsify our considerations. 

Finally, the resulting diluted DM abundances for the two scenarios HOT and COLD 
are plotted in Fig. [9j for the two different example values of M\ = 3 keV (~ lower bound 
on M l for e 1 = 0.001) and M l = 5 keV (> lower bound on M x for e x = 0.001). As 
indicated in the upper two panels, for a keVin with a mass close to the lower bound, 
one would need to correct the natural Xi abundance by at least a factor of about 35, 
which is already hard to achieve within the range of the hard BBN bound. However, it 
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Figure 8: Freeze-out temperatures and entropy production for the HOT and COLD sce- 
narios. The colour code is the same as in Figs. [5] and |6j 



is no problem if one only requires the weak BBN bound to be fulfilled. The lesson to 
learn is that the exact value of the bound is decisive for our framework, and it would be 
worthwhile to perform a detailed dedicated study. Requiring only the weak bound from 
now on, we can see from the plots that there is, especially for the COLD scenario, a sizable 
region in the parameter space, for which we can achieve the correct DM abundance. In 
particular, there is a certain trade-off between the parameters M l5 M 2 , and e^- Note that 
Mi can in principle be freely chosen, as long as the lower bound from structure formation 
is not violated, cf. Fig. |3j However, increasing M\ also increases the amount of entropy 
dilution, <S req , needed to obtain the correct DM abundance - which will cut the amount of 
the M 2 -62 space that leads to the correct abundance. In particular for the HOT scenario 
with Mi = 3 keV (upper left panel of Fig. [U]), we can see that there is a strong tendency 
to actually produce too much entropy for hot freeze-out of the Xi- in addition, this plot 
also shows that it would be worthwhile to study the intermediate freeze-out region of 
\2i which cannot be done in a semi-analytical way. Still, this could very well be the 
ideally suited parameter region to obtain the correct abundance. We want to stress once 
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Figure 9: Diluted DM abundances for the HOT and COLD scenarios. The colour code is 
the same as in Fig. [71 Note that the terms "HOT" and "COLD" refer to the freeze-out 
of X2- The different curves correspond to those shown in Fig. [HI 



more, however, that the true conclusion about this parameter region would require a more 
advanced numerical investigation of the Boltzmann equation. 

One bonus of our general treatment is that we are able to investigate a large class of 
models within our framework. By comparing the interaction Lagrangians with the Z- 
boson with the general forms in Eqs. ([1]) and (j2J), we can express our general parameters 
in terms of model paramaters by the following mapping: 



(ei, e2, 5, Mi, M 2 ) — > /(model parameters), (21) 

where / is a function that is characteristic for the concrete model under consideration. 
Hence, if we have a fully working set of parameters (ex, e 2 , 5, Mx,M 2 ), we can use Eq. (l2"Tp 
to obtain conditions on the model parameters. We have collected four such benchmark sets 
of parameters in Tab. [lj All these points lead to DM abundances that are in agreement 
with the 3a range for il DM h 2 obtained by WMAP-7 0. 

Note that, actually, it would be possible to generalize our considerations to models 
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Table 1: Four benchmark points that lead to the correct abundance and are all consistent 
with the Ly-a and the weak BBN bound. The classification "HOT/COLD" refers to 
the freeze-out of \2i and the mass M 2 is chosen such that the correct DM abundance is 
obtained. 

involving Z'-bosons, since this would in the zeroth order approximation only amount 
to a rescaling of the parameters ei, e 2 , and 5. However, in addition to that, it may 
be that the \% coupling to the Z-boson is not of pure axial vector form, cf. Eqs. (JTJ 
and (J2J), and/or that the function S, cf. Eq. 01]), is modified, together with the replacement 
(M Z ,T Z ) —> (Mz',Tz') in Eqs. (jl]), (J5J), and ([7]). However, in situations where one is not 
interested in the region near the resonance peak, and where an 0(1) approximation for E 
is sufficient, one can still absorb most of the new model aspects in a rescaling of (ei, e 2 , 5). 

In the following section, we will present some realistic cases where one could make use 
of our benchmark points. 



3 Example Models 

After having discussed the general framework of keVins, we now want to look at concrete 
examples to illustrate in which types of models a keVin could arise as DM candidate. We 
discuss two realizations of our framework: first we will give an example of keVins within 
the -G^SSM, with the Z-boson being the dominant mediator of the Xi annihilation. In 
that setup, the decay would involve a "flavour" changing coupling \2 with xi, and would 
hence be naturally suppressed compared to the annihilation channels, thereby justifying 
the use of our parameter 5. Secondly, we will illustrate how to extend our considerations 
the case of keV sterile neutrinos. This extension needs additional constraints to be taken 
into account, which impose new difficulties not present in our general setting. However, 
when carefully adjusting our parameters, also this type of DM could be accommodated. 



3.1 A SUSY example: £ 6 SSM 

The situation of DM in the £ 6 SSM ESI ESI ED] (see Ref. [EQ for a pedagogical in- 
troduction) is very similar to our general case presented in Sec. [2J we have one light 
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mass eigenstate singlino Xi, which could play the role of a WDM particle, and a heavier 
brother of that particle, the next-to-lightest singlino x 2 ; which can decay into light parti- 
cles. Both, Xi an d X2, are thermally produced in the early Universe, and the decay of X2 
is responsible for the entropy production. 

Our basic framework could actually be any supersymmetric version of the Standard 
Model, which will be denoted by SSM, as opposed to SM. Within this SSM framework, 
we need three ingredients: 

1. Two (in general distinct) SU(2)l doublet superfields with hypercharges +1/2 
or —1/2. As the notation suggests, these might be the two Higgs doublet superfields 
of the minimal SSM (MSSM), or of any extended model. 

2. Two distinct SU(2)l singlet superfields S1.2 with vanishing hypercharge. A generic 
example of such fields would be the singlets in the i? 6 SSM, which are usually referred 
to by a very similar notation as our singlets. 

3. Some quantum number that distinguishes the fermionic components si^ of Si^ from 
right-handed (SM singlet) neutrinos. This job would usually be done by i?-parity, 
which would in addition make the lightest inert mass eigenstate stable, although 
there could be other choicest 

In general, such a framework will tend to mix these gauge eigenstates (after electroweak 
symmetry breaking), which in particular applies to their fermionic components hi and 
Sj. We assume this mixing to be such that, at tree- level, hi mixes with Sj only by a 
small admixture of size while the mixing with Sj-n could even be further suppressed, 
paramaterized by an additional factor of S 1 ^ 2 . In such a situation, the mass eigenstates 
Xi,2,3,4 are given b}0 



( X1 \ 


/ 


1 


e l e /2 \ 


(s 1 \ 






-ei 


1 e /2 


hi 


X2 











w 


V 


_£l/2 


-e 2 1 J 


w 



(22) 



One should keep in mind that we are talking about a situation where the couplings to the 
Z-boson are suppressed, i.e., we have in mind that 6j <C 1. Such a mixing pattern may 

9 Note that this condition prevents the from mixing with neutrinos, which may by a big advantage 
or disadvantage of our setting, depending on the viewpoint: indeed, in the view of a model builder one 
avoids the strong constraint from the non-observation of the 7-line, while in the view of an observational 
astronomer one loses a striking signal. 

10 Note the order of the states on the left-hand side. 
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look somehow engineered at first sight, but it might very well be motivated in a concrete 
model: in the _E 6 SSM, this structure suggests itself due to the third generation singlet 
S3 being the only one to get a VEV. Note further that, since the matrix in Eq. (1221) 
mixes SU{2)l singlets with doublets, it is natural for the corresponding couplings to 
be proportional to the electroweak VEV, q oc v, which is normally small compared to 
any other VEVs. Furthermore, the "flavour-changing" vertices (i.e., Z~Xi~~X2 couplings) 
typically arise only at 1-loop level in the l^SSM, which suggests an additional suppression 
parametrized by our parameter 5 (or by £ = e^S). 

Due to their singlino nature, both x° have suppressed couplings to the Z. They are, 
however, charged under the full gauge group and couple with gauge strength to a heavier 
Z 1 gauge boson. Depending on the mass of the Z', these couplings could be suppressed. 
The best current lower limit on the mass of the Z'-boson has been obtained by the ATLAS 
collaboration [62], and it is given by Mz> > 1.83 TeV. On the other hand, there is no a 
priori reason in the -B^SSM for the Z' not to be even heavier, in which case the effective 
couplings of the Xi to the Z' are again very small. It is this limit that we will consider. 

The decisive step is to compare the couplings in the model to the general forms given 
in Eqs. ([I]) and ([2]). Starting with the flavour diagonal couplings, one obtains [58] : 

9 e i = -rr-Rzih where z = l,2. (23) 
iv 

Here, M z = \/ g 2 + g' 2 v is the usual Z-boson mass, Rza = 2 (m ± ) 2 yf cos2 @ ~ fi fi) 

is the coupling parameter to the Z-boson, and m + = ^tf are the "inert" chargino masses. 
Furthermore, s is the VEV of the third generation singlet scalar component S3, tan/3 
is the usual MSSM-like doublet VEV ratio, and the couplings fi, fa, and Aj arise from 
trilinear couplings of one singlino superfield to two Higgs doublet superfields in the ^SSM 
superpotential [52] . 

Using Eq. ( |23]) . one can express the general coupling ^ in terms of model parameters: 

W? cos 2 p - Of sin 2 (3 -, (24) 



where Ui = fi/Xi and cDj = fi/Xi- Taking the ratio between ei and 62 leads to 



11 = , 1~^ tan ^ » 1, (25) 
e 2 V ul - u\ tan 2 (3 v ' 

as fulfilled for all of our four scenarios. Hence, although the VEV ratio - determines the 
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scale of the couplings e«, one still needs some way to make the ratio in Eq. (1251) large 
enough. We will discuss two possible solutions below. 

Similarly, one can express the masses in terms of model parameters, 

Mi = i (26) 

2 

Like before, the size of the masses is determined by — , which suggests that both are 
significantly smaller than the electroweak VEV v. Taking again the mass ratio, 

— - = / 1 . 1 < 1. 27 

M 2 A 2 lo 2 uj 2 

In order to simultaneously fulfill Eqs. ( 125|) and ( |27|) . we have different possibilities: 

• flipped hierarchies: fi > f 2 , f 2 3> fi 

If Ai ~ A2, then this condition leads immediately to Mi <C M2, due to /1 <C j%. In 
turn, Eq. ([2"5]) simplifies to 

- - 1 = (28) 
• I ft -ft tan 2 

whose numerator is much larger than the denominator. 

• tuned tan/3: tan/3 ~ /2//2 

This assumption leads to u 2 — tan 2 /3 w in Eq. (pH|) . and hence to a very small 
coupling e 2 . However, this condition still allows all other combinations of couplings 
to be sizable. Thus, in order to also fulfill the condition in Eq. ( 1271) . it is necessary 
to impose an additional assumption that leads to fifi/Xi 72/2/^2- 

We still have to test our original assumption of a Z'-boson that us practically decou- 
pled: M 2 /v ~ 0.01-0.1, together with Eq. fl26]), suggests that s ~ 0(10-100) v. Hence, 
a typical annihilation cross section involving a Z'-boson instead of a Z-boson would be 
suppressed by a factor of about (M^/M^/) 4 , which is approximately (v/s) 4 < 10~ 4 . This 
has to be compared to the natural ej suppressions of the cross sections involving the or- 
dinary Z-boson, in order to ensure that the Z'-contribution is indeed subdominant. For 
e 2 ~ 10~ 4 , this leads to the condition 

M z , ~ s > - > O (IO 3 - 4 ) v. (29) 
£2 

But this and Eq. fl2ED would lead to M 2 < - v < 10~ 3 t; ~ 0.1 GeV, which is too 
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2.4- 1(T 2 
1.4- 10^ 5 


6.1 • 10" 3 
1.3 • 10" 5 




2.0 • 10" 3 
1.3 • 10^ 5 


h 

h 


3.5- 10- 1 
6.8- 10- 1 


1.0 
0.67 




1.0 
0.67 


Ai 
A 2 


1.1 
1.4- 10 _1 


1.0 • io- 1 

0.43 




8.4- 10" 2 
0.23 


5 

tan/3 
s [TeV] 


1.8 • 1(T 2 
0.52 
3.78 


1.8 • 10" 2 
1.5 
6.37 




5.3- 10~ 5 
1.5 
2.38 



Table 2: Parameter values for the i^SSM with approximate decoupling of the Z'-boson. 
Note that we have not been able to find a fitting parameter combination for Ce- 



small to be consistent with the reheating bound, cf. Sec. |2j A simple way out would 
be to assume that the Z'-boson mass receives additional contributions from further new 
physics, thereby invalidating the above assumption Mz' ~ s. Then we could have a much 
stronger suppression of the Z'-mediated couplings, justifying the small values for the 
parameter 5 from Tab. [IJ cf. Eq. j2]). This off-diagonal coupling arises in the E^SSM only 
at loop level, as it vanishes at tree-level for X a s ^> f a v, f a v [63]. Even though we cannot 
give an analytical formula, one could generically expect a size of 

Rzi2 ~ V Rzi\Rz22~^ (30) 

where 5 is some suppression factor. Comparing this to Eq. (|2j) leads to 

M z 

which translates into 5 = tt— S. Hence, 5 should be tuned to have a somehow similar 

2c w ' 

value as 5 in the benchmark scenario under consideration. The natural suppression factor 
for a loop would be about 5 ~ l/(167r 2 ), which is okay with scenarios Hd,c from Tab. [H 
while Ce,d would require some more suppression. That could come from the structure of 
the loop which could exhibit, e.g., a classical (61] or extended [63] GIM mechanism. 

Trying to find explicit realizations for the different benchmark points from Tab. [TJ it 
is easiest to perform a ^-minimization. Good fits were found for scenarios Hd,c and Cd, 
and they are listed in Tab. |2j Note that the fits for He and Cd both correspond to the 
tuned tan (3 case mentioned above, while the one for H D corresponds to the case of flipped 
hierarchies. 
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3.2 keV sterile neutrinos in a Left-Right symmetric context 

In this section, we would like to shortly illustrate that it is, in principle, possible to 
extend our considerations to the case of keV sterile neutrino Dark Matter, where the 
lightest sterile neutrino Ni is assumed to have a keV-scale mass. However, this extension 
also introduces additional constraints have to be taken into account, in particular the 
astrophysical bound on the radiative decay Ni —> wy [661 EZ] arid constraints from active- 
sterile neutrino mixing. The idea of using additional entropy production to dilute the 
thermal abundance of keV sterile neutrinos was, to our knowledge, first mentioned in 
Ref. [26]. The key point is to assume that the N\, although being a SM singlet, shares 
some interaction at high energies. A very concrete example of such a case has been worked 
out in Ref. [19], where a Left-Right (LR) gauge symmetry was assumed under which the 
right-handed (SM-singlet) neutrinos are charged. We will refer here to the notation used 
in Ref. [19J to give a flavour of what would change when extending our considerations. 

To generalize our considerations to the case of sterile neutrinos in an Li?-symmetric 
framework, we first make the identification (xi,X2) {Ni,N 2 ). Then, one immediately 
observes that the corresponding sterile neutrino masses M\ and M 2 are free parameters, 
and can be tuned to have, e.g., the values of one of our benchmark scenarios in Tab. [Til 11 ! 
Reproducing the correct values of the paramaters (ei,e 2 ,<5) is a bit less trivial. Still, 
it is possible to redefine our paramaters in terms of model paramaters: although the 
corresponding decay and annihilation processes are mediated partially by heavy charged 
Wft bosons, the structure of the diagrams does not change considerably. Because of the 
pure gauge interaction, we need €\ — e 2 , due to the gauge interaction, which is not possible 
in any of our scenarios. However, it is easy to find other points in parameter space where 
the correct abundance is reproduced. One such point with e\ = e 2 would be the following: 

(e l7 e 2 , 6, Mi, M 2 ) = (1(T 4 , HT 4 , 1, 2.6 keV, 4.4 GeV), (32) 

which produces an abundance of fl^h 2 = 0.113 and a reheating temperature that is in 
accordance with the weak BBN bound. In addition, the keVin mass of Mi = 2.6 keV is 
close to but still above the Ly-a bound. 

Taking this new benchmark point, we can compare the suppressed coupling to weak 
interaction from Eq. ([1]) to the strength of the suppressed coupling of Ni to the heavy right- 

11 Note that this statement is only true as long as there is no reason for the masses to have certain values 
in a particular model. There is still not too much work done in that direction, but among the possibilities 
to explain such a mass pattern involving a strong splitting are the split seesaw mechanism [33] (which 
can be supplemented by an An symmetry [34]), by L e — — L T symmetry [35] [36| (or other flavour 
symmetries |37j). by the Froggatt-Nielsen mechanism [38 s 39 , 40], or by the so-called extended seesaw [41] . 
as mentioned in Sec. [T] 
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handed HA^-bosons used in Ref. [19]. By comparing the cross sections, and disregarding 
the differences between the 75 and the Pr coupling for simplicity, we obtain 



e 12 = 0.0001 « 1.27 y/G F c w 



M W M Z 
Mw R ' 



(33) 



leading to Mw R = 2.97 • 10 5 GeV. This scale looks much smaller than the value of Vr ~ 
10 5 TeV quoted in Ref. [19], so apparently the numerical example given there does not 
correspond our benchmark point. This is, however, not too much of a wonder since we 
did not attempt to find all working parameter combinations, but rather wanted to give 
a proof-of-principle. We can go one step further by comparing our approximate decay 
width, Eq. ffT9]) . with the estimate Tn 2 > ^ f 9 2 [19] and again disregarding the 7- 
structure, which leads to a second generation active-sterile mixing angle of 9 2 ~ 6.9 ■ 10 -8 , 
in good accordance with the value obtained in Ref. pj)] . Note, however, that their scenario 
is somehow less restricted than ours, since the authors assume the gauge interactions 
to dominate the DM production, while they take the Yukawa couplings to be mainly 
responsible for the decay of A^. This disentangles certain requirements, which in our case 
can only be partially resembled by choosing a small value for 8. 

However, as already mentioned, additional constraints come into the game, the strongest 
being in this case the upper bound on the radiative decay width, r(iV 1 — > w~f), which arises 
from the non-observation of the corresponding monoenergetic X-ray line [66] [67] 1^1 Ac- 
cording to Ref. [19] , this bound would constrain the mixing angle 9\ to be less than about 
3.8 • 10 -4 for M\ = 2.6 keV. On the other hand, active-sterile neutrino mixing constrains 
a certain combination of masses and mixing angles to be larger than the square root of 
the solar neutrino mass square difference, 



supposed that we have a type I seesaw situation. It is exactly this bound, Eq. (134|) . which 
cannot be fulfilled with the masses and mixing angles obtained above, in agreement with 
the findings in Ref. [19]. Hence, indeed we have seen that it is the additional constraints 
that generically arise in the framework of keV sterile neutrinos, which could destroy the 
validity of an otherwise working benchmark point. Thus one has to be careful when 
extending our findings to this case, since additional strong bounds might play a decisive 
role. However, the constraint from Eq. ( I34"j) can be avoided in a seesaw type II situation, 
which naturally arises in Li?-symmetric extensions of the SM. This possibility has also 

12 Note that additional signals, like e.g. enhanced dipole moments [68], might show up in a cosmological 
context. 



M x 9\ + M 2 9\ > V8 ■ 10 5 eV ~ 0.009 eV, 



(34) 
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been worked out in detail in Ref. [19] and has the feature of being much more model- 
dependent, due to all the detailed considerations to be undergone when attempting to 
properly take into account neutrino mixing. 

Let us end by noting that, as mentioned in Sec. 12.41 the amount of entropy produced 
could be increased by involving more than one generation of decaying fermions. In the 
case of a keV sterile neutrino, this possibility is particularly attractive, as there would 
normally be yet another heavier sterile neutrino which could also contribute to the 
entropy production. In particular in the context of Li?-symmetry this particle must exist, 
and one could use it to dilute the N\ abundance even further, if required. 

So, indeed, it is possible to extend our general consideration to the case of keV sterile 
neutrinos, but the situation is more complicated and leads to a strong loss of generality, 
the latter being the main reason for us not to focus on the special case of keV sterile 
neutrinos here. 

4 Conclusions 

In this paper, we have proposed a simple model for Warm Dark Matter (WDM) in which 
two fermions are added to the Standard Model: stable "keVins" ( keV inert fermions) xi 
which account for WDM and their unstable brothers, the "GeVins" ( GeV inert fermions) 
X21 both of which carry zero electric charge and lepton number, and are (approximately) 
"inert" in the sense that their only interactions are via suppressed couplings to the Z. 
We have considered scenarios in which stable keVins are thermally produced and their 
abundance is subsequently diluted by entropy production from the decays of the heavier 
unstable GeVins. 

We have investigated in detail the question under which conditions it is possible to 
arrive at the correct relic abundance with these particles, which is done in our case by 
thermal overproduction and subsequent dilution of the abundance by additional entropy 
production in the decay of the \2 GeVin, which is assumed to have a mass around the GeV 
scale in our benchmark scenarios. We have identified four example benchmark scenarios 
for which it is possible to achieve the correct abundance, depending on the nature of the 
freeze-out of xi- 

This mechanism could be implemented in a wide variety of models, including E§ 
inspired supersymmetric models in which a keVin candidate naturally appears, or models 
involving sterile neutrinos. For example, in the E$ models the keVins are mainly Standard 
Model singlet fermions with only very small active admixtures due to mass basis rotations 
allowing them to have strongly suppressed couplings to the Z-boson. 

Although we have been able to give a proof-of-principle in this work, we nevertheless 
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have applied simplifications on several places, in particular in what concerns the freeze-out 
processes. The natural subsequent step would be to perform a more elaborate numerical 
study of the Boltzmann equation in order to extend our considerations to the intermediate 
freeze-out region of X2, and by this to directly calculate the freeze-out temperatures. 
Further investigations should be done, e.g., by identifying further regions in parameter 
space that lead to the correct relic abundance, and by using them to map other concrete 
models onto our general framework. We hope that we have been able to perform the first 
step in this direction and that many future works will investigate the framework of keVins 
in great detail. 
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